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An examination of the properties of many-electron conduction through spin-degenerate systems 
can lead to situations where increasing the bias voltage applied to the system is predicted to decrease 
the current flowing through it, for the electrons of a particular spin. While this does not necessarily 
constitute negative differential conductance (NDC) per se, it is an example of negative differential 
conductance per spin (NDSC) which to our knowledge is discussed here for the first time. Within a 
many-body master equation approach which accounts for charging effects in the Coulomb Blockade 
regime, we show how this might occur. 



I. INTRODUCTION 

Nonequilibrium properties like electronic conduction 
in molecular systems must be treated within a many- 
body nonequilibrium theory, and the extensive body of 
effects such treatments produce has attracted many re- 
search efforts on both the theoretical and experimen- 
tal sides pii2ii^i^i^i^i2ii^ In particular, models of transport 
through molecules and quantum dots have been shown 
to describe many nontrivial phenomena, of which the 
Coulomb Blockade effect is a well known example]^ Such 
models have been shown to describe cases where an in- 
crease in the source-drain bias on a small device coupled 
to macroscopic leads actually results in a decrease of the 
current through it i-'^-^i-'^^i-^^i-'^^i-^^i-'^^ This nonlinear behav- 
ior is known as negative differential resistance or conduc- 
tance (NDR or NDC), and its explanation must lie in the 
shifted states of the system and the switching of electron 
populations between them, but the exact mechanism may 
differ between the various cases. 

Several such mechanisms, many of which are actu- 
ally single-electron effects, are worth noting: the reso- 
nant double-barrier tunneling junction familiar in doped 
semiconductor workj^^ where an increasing bias pushes 
a resonant conduction state into the conduction window 
and then out of the conduction band of one of the elec- 
trodes, resulting in NDRji^ the case in which the elec- 
trodes themselves have narrow resonant features in their 
density of states, like an atomic-scale STM tip or an atom 
weakly coupled to a larger electrode, where the bias shifts 
the conducting levels of the electrodes into and out of 
alignment with each other j^^i^^i^^ the Coulomb-Blockade 
case where the bias charges the system in a way that 
kicks a level out of the conduction windowji^^>2a the 
more general case where the biasing actually conforms 
the molecule or causes a change in the interaction with 
phonons,— >^>22i>^ resulting once again in fewer available 
conduction levels. 

When the situation is compHcated by the lifting of spin 
degeneracy, the interplay between the occupation of spin 
levels and their coupling to the leads can result in spin- 
dependent effects as well. This has often been explored in 
cases where the leads are ferromagnetic ^^^^^ for example 



in the spin-blockade or spin field effect transistor More 
recently, spintronics without polarized leads have been 
suggested,.28^^^^^^^'26^^^^ where the leads 
would generally contain electrons of both spins, which 
would scatter through the system with different trans- 
mission properties. In a practical application the device 
might perform various transformations on spins rather 
then just act as a current switch. Thus it makes sense to 
develop the concept of conduction or resistance per spin, 
with the understanding that the same range of nonlin- 
ear phenomenon that is of interest for the total current 
can occur here for spin-dependent current. Specifically, 
the study of NDC is naturally complemented by NDC per 
spin, which is exhibited whenever increasing the bias volt- 
age on a conducting device causes the current through it 
for one spin to decrease, while the total current does not 
necessarily decrease. 

In this paper we take an illustrative look at a novel 
mechanism for the phenomenon of negative differential 
spin conduction (NDSC). We show that in cases where 
the charging of a quantum dot is a dominant energy 
scale of the problem and the spin degeneracy is lifted, 
NDSC can occur. The basic mechanism involves popula- 
tion switching between the two spin levels. In Section ITTl 
we describe a simple device in which NDSC may appear 
and be of interest, and explain the multi-electron master 
equation approach we employ for calculating the spin- 
polarized current. In Section ITTTl we display and analyze 
the results. Finally, we discuss our conclusions in Section 

llYl 

Some of the topics touched upon in this work have 
also been addressed by Raphy Levine over the past 
decade.— It is a great honor to dedicate 
this work to him, on the occasion of his 70*'' birthday. 



II. MODEL 

Perhaps the simplest and most abstract spintronic 
device one might imagine consists of a system with 
a single (energetically relevant) electron level coupled 
to two metallic leads. By making this level non- 
degenerate in the spin degree of freedom in any de- 
sired way, one can achieve filtering behavior by tun- 
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ing the conduction window so as to contain only one 
of the spin levels.— i^>^i^>^ If the system is small, 
one expects that the charging should become an impor- 
tant energy scale in the problem, and the population of 
electrons the device contains at any given time should 
not vary greatly from the neutral number of electrons. 
At this regime a many-electron master equation treat- 
ment^i2id8d9i50i5i|52^^^ can be expected to provide 
a good approximation of the dynamics, particularly for 
larger voltages. 

We consider a model for a spin-filter device which is de- 
scribed by a two single spin levels. Such a device can be 
an atom, a molecule, a quantum dot, or any system with 
discrete levels that are well separated. The single elec- 
tron levels on the device are coupled to two leads with 
chemical potentials /zl and and coupling constants 
lL(R),ij [i — o-jb), which we will assume to be equal to 
lL{B),ij — jSij. The model Hamiltonian of such a de- 
vice (not including the leads, since the current and level 
populations will be calculated within a standard multi- 
electron master equation approach}^ is given by: 

HD = {e + 5) a^a + [e - 5) b^'b + A {a^a + b^'b - Nof , 

(1) 

where a and b are single particle annihilation operators 
corresponding to the two spin levels (a =1 and b =1, 
respectively) and A^o is the neutral number of electrons. 
We also introduce the spinless level energy e, the spin 
energy shift S and the charging energy A. We note in 

I 



{tpi\G\ipj) = gij 



Multi-electron effects will be considered only in the form 
of charging energy. Since these values will be used in a 
rate-process calculation rather than a full quantum for- 
mulation, constructing the multi-electron states them- 
selves is actually redundant, and Eqs. (|3bp-(|3dp along 
with the single particle data will provide all the neces- 
sary information. 

The transfer rates between the multi-electron states 
are given bys^ 

Ri,a^l3 = ^'^^ Qii3, (4) 

where the four multi-electronic states are labeled by the 
Greek indices |a) = |'T-a"^?i["^) and = \ni^^nf ^), such 
that (for instance) |00) is the empty state, |01) is the 



passing that the notation is only for convenience and no 
real assumptions are made as to the symmetry of the 
shift. In fact, the results are relevant to any two sep- 
arate channels with different energies, for instance two 
quantum dots of slightly different energies each of which 
is coupled to different leads with a charging interaction 
between them. 

We now describe the approach taken to construct the 
multi-electron master equation, suitable for the above 
model, from single electron data. If one neglects spin- 
dependent multi-electron effects, then it is formally 
straightforward to build from a set of one-electron Hamil- 
tonian and spin eigenfunctions an anti-symmetric basis 
of multi-electron wavefunctions. Limiting the discussion 
to only two levels, one can define: 

^-Kina = n "Pn,- (2) 

Here A12 is the two particle anti-symmetrization operator 
and the states are identified by their (spin-dependent) 
level occupations Ui (0 or 1 for fermions) . Using this anti- 
symmetric multi-electron wavefunction we can uniquely 
and conveniently determine the nonzero matrix elements 
of a general many-body operator G required to construct 
the master equation. According to the Slater-Condon 
rules where only single electron integrals are taken into 
account: 

(3a) 

(3b) 

(3c) 
(3d) 
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state where only level b is occupied, and |10) is the state 
where only level a is occupied. The lead index in the 
above is ^ e {L, R}. For reasons that will become clear 
below, we also define the total transfer rate summed over 
both leads: 

Ra->l3 = ^ Re,a^l3- (5) 
i 

Following the Slater-Condon rules (cf., Eqs. l(3c|) and 
(|3d|) ). the coupling between the multi-electron states, 
is related to the single electron level coupling (or 
the imaginary part of the self-energy)^ and is given by 
^e,af3 = = CL,b) if the two multi-electronic states 

differ only by the occupation of level i, Vt^ap = 
if they differ only by rii and rij and rii — rij = 1, and 
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i=i 

= 5l2<5n2-r4-l'^rii-ni+l + 521<5n2 -n^ +l<5ni -n'^ - 1 • 



r^,a/3 = otherwise. As noted above, 7^^^ is the ma- 
trix element of the single electron level coupling. in 
Eq. |[4)) is related to the Fermi-Dirac function, /(e): 

= ^ 1 - /(fa - e/3 - < Np, (6) 

where Na — J2i ""-i^^ is the number of electrons in state 
a. The state energies are calculated from the Hamil- 
tonian ^ and amount to AN^, {e + 6) + A {1 - Nof , 
(e - (5) + A (1 - iVo)^ and 2e + A (2 - Nof respectively 
for the states |00), |10), |01) and |11). 

Once the rates are known, the linear master equation 
system can be read from the detailed-balance condition 
for steady-state: 

Ra^pPa - ^ Rp^aPp = 0, (7) 

where Pa is the probability that the system is in a multi- 
electron state a. The current at steady state is given 
in terms of the steady state occupation probabilities and 
can be expressed as^,; 



h = —e'Y^ Re,a^pPaSal3, 
af3 



(8) 



where 



Sa/3 



+ 1 Na < Np, 
-1 Na > Np, 
Na = Np. 



(9) 



Intuitively, this expression states that current flows out 
of lead £ whenever an electron flows from it into the 
device, with the inverse also true. Following a similar 
line of physical reasoning leads to an expression for spin- 
polarized current: up or down current flows out of lead 
£ whenever an up or down electron flows from it into 
the device. Assuming no coupling between levels with 
different spin, the spin-dependent current is given by: 



and 



Sa(f))a/3 



h,a(b) = — e Ri^a^pPaSa(b)ap, 
a/3 



-1(0) Sa<SpANa<Np, 

0(+l) Sa> SpANa<Np, 

0(-l) Sa< SpANa> Np, 

-1(0) Sa> SpANa> Np. 



(10) 



(11) 



Here, Sa = J2i ■'i"^ where si"' — ±1 for spin up (a) or 
down (6), respectively. 

The linear master equations can be solved analytically, 
but the form of the solution is rather cumbersome and 
has no real beneflt. They can also, of course, be solved 
numerically, which is the method chosen for this work. 
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Figure 1: The current (upper panel) and state populations 
(lower panel) as a function of bias voltage, for the parameters 
P = 20/7, fi — 0, e = 'y, S — 7/5 and A = 7/2, where the 
neutral number of conduction electrons has been taken to be 
A''o = 1 and the voltage is applied symmetrically. Note the 
distinct region of negative differential conduction (decreasing 
current) for each spin current, which always occurs simulta- 
neously with a decrease in population for a state populated 
with the matching spin. 



III. RESULTS AND ANALYSIS 



As the above model is more or less identical to the one 
commonly used to explain the Coulomb blockade, it is 
no surprise that plotting the I — V characteristics of the 
system shown in the upper panel of Fig. [T] immediately 
displays the well-known nonlinear current steps typical of 
Coulomb blockade effect.^" The model parameters used 
here are (3 — 2O/7, = 0, e = 7, 5 = 7/5 and A — 
7/2. The neutral number of conduction electrons has 
been taken to be TVq = 1- Conduction peaks or rises 
in the total current are expected in this formalism when 
there exists an energy difference between two states with 
electronic occupations that differ by one {Na — Np = ±1), 
which is also the energy of an electron occupied in one 
lead but not the other (when AE = Ea — Ep = ^±eV/2). 
In other words, the total current can rise at any bias 
voltage Vb where the conduction window is expanding 
so as to contain some spectral line of the system. For the 
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present model, this occurs when 

\e±5 + ^[(2-Naf -{l-Naf). 

(12) 

This is clearly the case for the current shown in the upper 
panel of Fig. [H where the four steps observed in the total 
current appear at O.67, I.47, 2.67, and 3.47 correspond- 
ing to transitions between states |00) ^ |01), |00) ^ |10), 
|01) <-> |11), and |10) <-> |11) , respectively and according 
to Eq. (fT2|) . 

Turning now to discuss the current per spin also shown 
in the upper panel of Fig.dJ we still observe the Coulomb 
blockade steps, however, the direction of the step can be 
either positive or negative. This is an example of a neg- 
ative differential spin conduction where an increase of 
the bias voltage is followed by a decrease in the current 
per spin. The NDSC occurs for both spins in this case, 
and at a different bias voltage for each spin-current. The 
first drop in current occurring for spin type b is also ac- 
companied by a sudden drop in the population of the 
|01) state (shown in the lower panel of Fig. [TJ, as the 
change in chemical potentials begins to allow the pop- 
ulation of the 1 10) state. This switching of populations 
between the states is reminiscent of another example of 
nonmonotonic changes in occupation predicted to occur 
in a system of two electrostatically coupled single-level 
quantum dotsjS"^. 

There is a simple "hand waving" explanation for such 
behavior: the current for each spin consists at low bias 
of contributions proportional to the probability that the 
system is in some state a and to the rate of transitions 
between state a and state |00) (in which none of the 
states are occupied), where a can be either |10) for spin 
up current or |01) for spin down current. We therefore 
expect that at any chemical potentials where the relevant 
Fermi functions and hence the rates are nearly constant 
at the relevant energy, the current will, to a good approx- 
imation, be linearly proportional to the population. The 
population, in turn, decreases whenever the shifted ener- 
getics allow the occupation of a new state. For this to be 
possible, the bias must be applied in such a way that not 
all states become occupied simultaneously. This is the 
reason the central chemical potential has been placed 
below the levels. The second current drop in the figure, 
which occurs at higher voltage and for the a spin cur- 
rent, can be explained by a similar argument - this time, 
however, the depopulation of the |01) state is the one 
involved. 

It is worth pointing out that the population shifts are 
such that in regions where the chemical potentials are 
far from any levels, any states that energetically can be 
populated become so with equal probability. At low bias 
only one state has the entire population, then as the bias 
is increased the population is shared equally between two 
states, then between three, and finally between all four 
states. This is the cause of the downward shifts in the 
population which result in the NDSC. For the example 



shown in Fig. [T] conduction sets in when the population 
of state 1 00) decreases from its maximal value of 1 until 
both states |00) and |01) are equally populated. Then 
NDSC occurs when both states |00) and |01) lose popu- 
lation to state 1 10) until all three state become equally 
populated. It is also evident that in systems with more 
electronic levels, NDSC due to population switching will 
become weaker if the separation between the states is 
small. Noticing this fact also clarifies the role of charg- 
ing in NDSC, as without charging all the states which 
include the same energetically occupiable levels would 
become populated at the same bias voltage. 

The theoretical phenomenon of NDSC and its physics 
are easy to understand, and the mechanism we suggest 
for it here simple, but two important questions remain: 
when will it occur, and how can it be observed experimen- 
tally? Answering the first question formally is a matter 
for the analysis of the expressions for the spin currents: 
by taking their derivatives and looking for a local maxi- 
mum in the voltage, an exact condition could be worked 
out in principle. In practice, the analytical development 
involves the solution of nonlinear equations, may or may 
not be possible and of interest, and is beyond the scope of 
the present work. Instead, a look at a part of the surface 
of transition in parameter space between regions where 
NDSC does and does not appear (see Fig. [2]) is enough to 
convince oneself that the exact conditions for NDSC are 
nontrivial. If one is more interested in the approximate 
limits where the drop in the current constitutes a siz- 
able fraction of its maximum value and where the master 
equation is vaHd, these can be expected when the tem- 
perature is smaller than the splitting between the spin 
states, > 1; when the conduction resonances are nar- 
row enough such that ^,5 > 7; and when the charging 
energy is of the order of the level spacing, A e > 7. All 
criteria can be met, for example, for systems of nanome- 
ter dimensions where the charging energy and the level 
spacing can be tuned by simple changing the size. Also, 
as mentioned above, some asymmetry in the application 
of the chemical potential and/or bias voltage is required, 
as the spin dependent effects happen to cancel out com- 
pletely when the chemical potential is exactly between 
the energies of the two levels and the bias is applied 
symmetrically. In addition, NDSC requires some state to 
become energetically occupiable at a bias voltage higher 
than one at which a spin current exists. More precise 
conclusions require a calculation similar to the one done 
to produce Fig. [21 which takes negligible computational 
effort and can be easily extended to more detailed sce- 
narios. However, the effect clearly occurs for an extraor- 
dinarily wide range of parameters, as can be seen in the 
figure. 

Addressing the second question posed above, pertain- 
ing to experimental observability, the following is pro- 
posed: NDSC is obviously equivalent to NDC whenever 
the spin components of the current are observed sepa- 
rately. This can be achieved in any experiment where in 
addition to fiowing through the system described here. 
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Figure 2: A plot of a part of the surface of transition between 
regions where NDSC does and does not occur for spin b is 
shown, for e = 7 and /i = 0, with the voltage applied sym- 
metrically. Note that at the limits of high temperature, small 
charging energy and small level splitting there is no NDSC, 
but that in general the behavior is complex. 

the current also flows (while retaining spin coherence) 
through some sort of spin beam splitting device which 
separates it into spin components before the current is 
measured. Such devices have been suggested in previous 
works . ^^1^° 



model of a quantum dot coupled to metallic leads, we 
have pointed out a mechanism that gives rise to negative 
differential spin conduction. The fact that NDSC occurs 
in such a basic model for a wide range of parameters sug- 
gests that it represents a real physical phenomenon. We 
have also discussed when effects of this type can be ex- 
pected to occur (the temperature must be low enough, 
conduction peaks narrow, and charging should be sig- 
nificant as this is strictly a many-particle effect), and 
suggested how an experiment in which they might be 
measured could be carried out. 

The NDSC effect, as caused by population switching or 
any other mechanism, is similar to and in special circum- 
stances identical to the NDC effect which has been ob- 
served in a variety of nano- and meso-scale experiments. 
It is characterized by an increase in bias voltage over a 
junction resulting in the decrease of the current for elec- 
trons of one particular spin. One way to observe NDSC 
directly is to measure the current after it passes through 
a beam splitter. Regarding the population switching 
mechanism: if the energetics are tuned so that increas- 
ing the bias allows the sequential occupation of several 
states, which with charging limiting the total occupation 
results in a nonmonotonic behavior of the populations, 
NDSC can be caused by a decrease in the population of 
a state which is instrumental in the conduction of one 
spin. When this happens without significantly effecting 
the electronic fiow rates to and from that state, a de- 
crease occurs in the contribution to the current from the 
term proportional to the population, leading to NDSC. 
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